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Previously  the  Legendre-Hermite  polynomial  expansion  of  Boltzmann's 
equation  (1)  has  been  used  to  calculate  the  velocity-field  charac¬ 
teristic  of  InP  at  room  temperature  (2).  In  this  report  application 
of  the  method  to  a  millimetre  wave  InP  Gunn  diode  is  described. 
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An  InP  n  -n-n  plane  diode- with  a  base  width  of  0.7  um  and  doping  of 
16  3 

3x10  /cm  has  been  studied  at  an  assumed  operating  temperature  of 
493  K.  The  exact  doping  profile  is  included  in  Fig.  3. 

These  data  were  chosen  because  a  similar  diode  was  simulated  at 
the  University  of  Lille  by  Friscourt  et  al.  (3)  using  an  energy 
relaxation  method  and  it  was  thought  interesting  to  compare  both 
methods. 

In  our  simulation  method  carrier  transport  was  calculated  by  the 

LHE  method  and  Poisson's  equation  was  solved  simultaneously,  taking 
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time  steps  of  about  10  sec.  The  boundary  conditions  at  the 
contacts  were  cyclical,  that  is,  electrons  leaving  the  device  through 
one  contact  were  injected  with  the  same  velocity  at  the  other.  To 
make  sure  the  electrons  in  the  anode  contact  layer  have  enough  time 
to  cool  off  both  contact  layers  were  taken  fairly  thick. 

Fig.  1  shows  the  d.c.  current-voltage  characteristic.  An  operating 
point  of  3  /  was  chosen  and  a  sinusoidal  a.c.  voltage  at  a  frequency 
of  143  GHz  was  superimposed.  The  resulting  current  waveform  was 
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Fourier  analyzed  and  the  power  output  at  the  fundamental  frequency 
calculated.  This  is  shown  in  Fig.  2  as  a  function  of  a.c.  voltage 
amplitude.  At  the  point  of  maximum  power  output,  1.5  V  a.c.  ampli¬ 
tude,  pictures  were  taken  of  the  carrier  density,  field  and  velocity 
profiles  at  eight  points  in  the  a.c.  period.  These  are  shown  in 
Fig.  3.  Clearly  the  device  is  operating  in  an  accumulation  layer 


transit  mode.  Also  it  can  be  noted  that  the  carrier  velocity  over¬ 
shoots  the  maximum  static  velocity  in  a  large  portion  of  the  device 
during  part  of  the  period. 

Finally  the  behaviour  of  the  diode  in  an  oscillator  circuit  was 
studied.  The  circuit  and  the  resulting  current  and  voltage  wave¬ 
forms  are  shown  in  Fig.  4.  After  an  initial  transient  caused  by  the 
abrupt  turn-on  of  the  bias  voltage  the  circuit  settles  into  a  stable 
oscillation  condition.  It  can  be  seen  that  the  device  current 
contains  an  appreciable  second  harmonic  component.  Also  the  average 
current  is  lower  than  in  the  non-oscillating  condition  (cf.  Fig.  1) 
which  suggest  application  as  a  self-oscillating  mixer. 
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Fig.  3a.  Electron  density  (cm  )  and  valley  population  (relative) 
as  a  function  of  position  (in  um)  during  an  a.c.  cycle.  The  text 
at  the  top  of  the  pictures  gives  the  phase  of  the  a.c.  voltage. 


V 


4b.  Device  current,  excluding  capacitive  current. 
Vertical  scale  100  mA/div. 
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